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ABSTRACT 



Aims. To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular 
.^r cores, whose asymmetries trace infall and expansion motions. 

O Methods. The high-mass star forming region SgrB2(M) was observed with HerschelfUlFl (HEXOS key project) in various lines of HCN and 
its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density 
(-H power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. 

f^ ^ Results. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity 
I (HCN to H"CN). This is most evident in the HCN 12-11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field 
O changes from infall in the outer part to expansion in the inner part. 

}--i Conclusions. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates 
C/3 _ in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally 
C^ , disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HlFl were critically important. 
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^' 1. Introduction crease due to ionization, stellar winds, bipolar outflows) push 

, ^.^ the gas outwards. Finally, this dominates over gravitational ac- 

^ Massive stars form deep inside molecular clouds. While this ob- ^^^^^^^^ ^^^j ^^e molecular cloud is disrupted. 

t^ scures the process at infrared or shorter wavelengths, it offers the j^^ ^-^^^ molecular cloud Sagittarius B2 is the most mas- 

O possibility of investigating high-mass star formation using dust ^j^^ ^„j ^^^^^^ star-forming region in our Galaxy. It is located 

O emission and molecular lines from the radio to the far-infrared, ^j^^^ ^^ ^^^ q^I^^^j^ ^^^^^^ ^^ ^ ^-^^^^^^ ^f 7 8 j^p^ (ReidHH 

T-H . The shapes of molecular lines are a particularly powerful tool, 2009). SgrB2(Mai n) is a t ight cluster of ultracompact Hii re- 

^ ; since they trace the density and temperature gradients as well as gj^^^ ( Gaume & ctliiiia il990 ). associated with a massive 

the velocity field. Even when the cores are not spatially resolved, j^^j molecular core. Its luminosity is estimated to be 6.3 x 10^ 

radiative transfer modehng of a variety of hnes can constrain l^ ^Goldsmithelali \m^, the most luminous hot core of the 

H , their structure. The unique spectral coverage of HIFI opens up Q^^^y ^^^^ j^f^n ^„^ ^ decelerating outflow were seen by 

. . . the full potential of this technique, since both the excitation InTn et al J d2008h 

of rotational transitions and the dust opacity increase with fre- ^^ ^^-^ ^^; ^^ ^^^^^^ ^^^ ^^^^j^^ ^^^ ^j^^ ^^^^^. ^^j^ 

quency. While the former allows us to study higher temperatures ^^ SgrB2(M), based on Herschel/niFl observations of HCN and 

and densities, the latter leads to stronger continuum radiation, complemented by APEX data, 
thereby facilitating the production of absorption features. 

A particularly important stage of high-mass star formation is 

the hot core phase. When massive stars begin burning hydrogen, 2 Observations 
they are still accreting and surrounded by dense molecular gas, 

which is heated and emits a rich spectrum of molecular lines. Herschel is a new space observatory with a 3.5 m antenna operat - 
Various feedback mechanisms (radiation pressure, pressure in- ing at submm to far-infrared wavelengths (Pil brattet am2010l) . 
Its spectral line receiver is HIFI (Heterodyne Instrument for the 



* Herschel is an ESA space observatory with science instruments Far Infrared, de Graauw et al. , 2010). The guaranteed time key 
provided by European-led Principal Investigator consortia and with im- project H EXOS (J/erac / 3e//HI FI Observations of EXtraOrdinary 
portant participation from NASA. Sources, iBergin et al.L l2010l) includes a full line survey of 
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Fig. 1. HerschelfHlFl observations of HCN in SgrB2(M). The green model has a constant infall, the red model has infall in the outer 
part and expansion in the inner region. Note the different asymmetries in the lines, which are reproduced by the red model. Upper 
state energies are stated above the plot. 



SgrB2(M) (ffyzooo = 17''47"'20.35' and dyzooo = -28°23'03.0"). 
In March 2010, this was started with HIFI bands la, lb, 2a, and 
4b, providing coverage of the frequency ranges 479 - 726 GHz 
and 1051 - 1121 GHz. Each band contains one rotational tran- 
sition from HCN, H'^CN, and HC'^N (Fig. [B, and two rota- 
tional transitions within their first vibrational state, which lies 
around 1000 K higher in energy (Fig.[2li. Trans ition frequencies 
of HCN (Thorwirthetal., 2003; Fuchs et al., 2004) are taken 
from CDMS (Mulleret al., 2001. 2005) . HIFI spectral scans are 
carried out in dual beam switch (DBS) mode, where the DBS 
reference beams lie approximately 6' apart. The wide band spec- 
trometer (WBS) is used as a back-end, providing a spectral res- 
olution of 1 . 1 MHz over a 4-GHz-wide intermediate frequency 
(IF) band. To allow separation of the two sidebands, bands la, 
lb, 2a, and 4b were scanned with a redu ndancy of 4, 8, 8, and 
4, respectively ("Comito & Schilk^ l2002h . The data were cali- 
brated using the standard pipeline released with version 2.9 of 
HIPE doTtlbOlOh . and subsequently exported to CLAS:fl using 
the HiClass task within HIPE. Deconvolution of the DSB data 
into single-sideband (SSB) format was performed using CLASS. 
All the HIFI data presented here, spectral features and contin- 
uum emission, are deconvolved SSB spectra. The horizontal and 
vertical polarizations differ by around 10%, and were averaged. 
The intensity scale is main-beam temperature, and results from 
applying a beam efficiency correction of about 0.68. 



' Continuum and Line Analysis Single-dish Software, distributed with 
the GILDAS software, see http://www.ira m.fr/IRAMFR/GILDASi 



APEX (Atacama Pathfinder Experiment) is a 12-m telescope 
located at 5100 m altitude in the Atacama Desert, Chile, one of 
the optimal sites for submm observations on Earth ('Gliste n et all 
2006). The lines presented here (Fig. [3]) are from a study of the 
structure of 12 hot cores (Rolffs et al. 2010, in prep.). Since the 
baseline levels in the ground-based data are not determined reli- 
ably enough to extract th e continuum, the 85 0yum continuum is 
taken from ATLASGAL dSchulleret al.L 120091) . 



3. Modeling 

The spherical Monte C arlo radiative transfer code RATRAN 
(Hogerheiide & van der Tak . 2000) was employed to compute 
the molecular lines and the dust continuum. Using an itera- 
tive by-eye comparison, the fit to the data was obtained in two 
steps. First, the physical structure (density and temperature) was 
adapted to match the continuum. Second, the HCN abundance 
and the velocity were optimized for the lines. 

The model we present here is centrally heated, its density 
following a radial power law with index 1.5. The dust opacity 
was taken from Ossenkopf & Henning ( 1994), without either ice 
mantles or coagulation. At the inner radius, 1 340 AU, the density 
is 1.8x10^ cm"-' and the temperature is 1500 K. The temperature 
was computed in an approximate way from central heating, us- 
ing the diffusion equation in the inner part and balance between 
heating and cooling in the outer part, assuming an effective tem- 
perature for the dust photosphere of 58 K. The temperature has 
a steeper gradient in the inner part, and falls off to 20 K at the 
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Fig. 2. HerschelfHlFI observations of vibrationally excited HCN 
in SgrB2(M), overlaid with the pure-infall (green) and the expan- 
sion+infall model (red). The kinematics are not well constrained 
in this hot, innermost region. 





'| M i|iiihi5 : '>l'i'l 



ftH'XN 4-3 
345.3 GH2 




M |ll H l !. l| M l| 



HC"N 4-3 
344,2 GHz 



'Ij'' y 




Fig. 3. APEX observations of HCN in SgrB2(M), overlaid with 
the pure-infall (green) and the expansion+infall model (red). 
Although less clear than in the HIFI data, the same asymmetry 
change can be seen. 
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Fig. 4. Sketch of the scenario described in the text. In the outer 
parts of the cloud, infall dominates. In the inner part, feedback 
has set in, in this example envisioned as multiple outflows, and 
gas is expelled from it. At the same time, infall continues along 
certain paths. 

Table 1. Continuum fluxes of SgrB2(M) and the model. 



Frequency 


Beam size 


SgrB2(M) 


model 


[GHz] 


["] 


[K] 


[K] 


345 


18.2 


3.8 


3.5 


500 


43.1 


2.4 


2.4 


600 


35.9 


3.8 


3.8 


700 


30.8 


5.7 


5.7 


1100 


19.6 


14.6 


14.2 



Note: The first row is the 850yum flux from LABOCA on APEX, which 
contains also lines. The next rows are the baseline levels at the given 
frequencies in HIFI bands la, lb, 2a and 4b. The baseline level of HIFI 
is a reliable continuum measurement, since it is as strong as saturated 
absorption lines and compares well with ground-based continuum data. 



exact shape of this, particularly the transition from infall to out- 
flow, is not constrained by the observations, hence modeling the 
true velocity field by adopting this spherical approximation is an 
oversimplification (see Fig.|4]i. 

A comparison with the observations (Figs.[T]|2l and[3]l shows 
that the red model reproduces more closely the observed fea- 
tures. In particular, to reproduce the changes in the asymmetry 
that occur both from lower to higher J and from higher to lower 
optical depth (from HCN 6-5 to H'^CN 6-5), a change in the 
velocity field is required, from infall in the colder, outer parts to 
expansion in the warmer, inner parts. 



outer radius of 5.6 pc, where the density is 7x10^ cm ^. The 
model reproduces the continuum very well (Table [1]). Its cen- 
tral dust optical depth is 0.4 at 345 GHz, 0.8 at 500 GHz, 1.5 at 
700 GHz, and 3.7 at 1.1 THz. 

To compute the lines, we used a turbulent 1/e half width of 
7 km/s (12 km/s FWHM). The lines were additionally broad- 
ened by their high optical depth. We found no evidence of sig- 
nificant variations in the turbulent width with radius. The HCN 
abundance increases with temperature from 10"*^ at temperatures 
below 100 K, to 3x10-** between 100 and 300 K, to 3x10"^ at 
temperatures above 300 K. The first jump is mainly needed to re- 
produce the intensities of the isotopologues and of vibrationally 
excited HCN, the second jump serves to fit vibrationally excited 
H'^CN. These jumps could be due to evaporation of ice mantles 
and/or increased chemical p roduction. We u s ed a ' ^C/'^C ratio 
of 20 and '^N/'^N of 600 (IWilson & Rood 1 1994 . The green 
model has a constant infall velocity of 2 km/s, and the red model 
has the same velocity in its outer parts, but an expansion in its 
inner parts. Figure |5] shows the velocity field. We note that the 



4. Discussion 

We now discuss the velocity field and how it is physically deter- 
mined. 

In a gravitationally contracting sphere with density power- 
law index p, a first guess and upper limit to the infall speed 

would correspond to free fall, i.e., Vfi(R) = J — ^^, where the 
mass inside the radius R is M{R) oc R^^p (for p < 3), such that 
Vjf ex /?T~. In our model, this free-fall speed is roughly (at a 
radius of 1 pc) 10 times higher than the 2 km/s of the model. 
Such high infall velocities can be excluded because they would 
produce lines that are far more asymmetric than observed. 

The mass accretion rate M - Vacc(^) j^ . so Vacc ^^ R''^ 
for a constant mass accretion rate. Hence, infall accelerates in a 
spherical steady-state model with p = 1.5. However, an accel- 
erating infall, at least in the inner part, is inconsistent with the 
data. For simplicity, we chose a constant infall velocity, corre- 
sponding to accretion rates onto the cluster of 3x10 ' Mo/yr at 
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Fig. 5. Physical structure of the model. Density is plotted as 
dashed green, temperature as dotted blue. The velocity field of 
the infall+expansion model is shown as the solid red curve (the 
infall model has a constant velocity of -2 km/s). At the transition 
between infall and expansion, the temperature is 65 K and the 
density is 6x10'' cm""*. 



the outer boundary and 10"' Mo/yr at a radius of 10^ AU. The 
excretion rate reaches 8x10"^ Mo/yr at 2x10^* AU. 

The real velocity field must be non-spherical, and infall 
could we l l con tinue in the inner part (see Fig. |4] and also 
lOin et all (l2008l) \ The peak on top of the HCN 12-1 1 absorp- 
tion (Fig.[l]) could be a hint of that infall. Since the H'^CN tran- 
sitions are less asymmetric than modeled points in that direction, 
one may suspect that a linear combination of the pure infall and 
expanding models would produce the line spectra very well. For 
the HCN lines, because of their high optical depths, such a linear 
combination of intensities does not naturally reproduce the ob- 
served spectra; more sophisticated modeling, in which we aban- 
don the spherical approximation, would be required. For over- 
lying model and data (Figs. [T] [2] andO, we assumed a source 
velocity of 64 km/s, although 60 km/s would reproduce more 
closely the high-J lines. This blue shift of absorption towards 
higher J is consistent with our picture of a velocity field chang- 
ing from infall to expansion, and is indicative of a deviation from 
spherical symmetry. 

While gravity is clearly responsible for the infall, there must 
also be feedback mechanisms from newly formed stars working 
against gravity, thus regulating star formation. The velocity of 
the infalhng envelope must be lower than that of pure gravita- 
tional collapse by around 90%. Apart from magnetic braking, 
turbulence induced by bipolar outflows may be such a long- 
range feedback m echanism providing the necessary pressure 
(IWang et al.Ll2010i) . The massive stars in the center dominate the 
luminosity, and drive the expansion by means of stellar winds, 
radiation pressure (including reprocessed light), and heating and 
ionization of the gas, in particular thermal pressure of their Hii 
regions (Krumholz & Matzner, 2009). Since in our model only 
1% of the total mass of 7x10^ Mq is expanding, the outward mo- 
mentum remains much lower than the inward momentum, indi- 
cating that feedback is not yet able to disrupt the cloud. 



5. Conclusions 

The changing asymmetry of the HCN lines clearly indicates 
a reversal of infall. Expansion motions dominate in the inner. 



warmer parts, and infall dominates in the outer, colder parts of 
the core. This can be naturally explained by the onset of feed- 
back from massive stars. Our radiative transfer modeling demon- 
strates the power of constraining the source structure by fitting 
line shapes. A model that successfully reproduces all the fea- 
tures cannot have spherical symmetry, and the exact geometry 
probably has to be constrained by interferometric observations, 
to determine the free parameters. Although the HerschelfHlFl 
data do not have the desired spatial resolution, they place strong 
constraints on any such model, since they have to be reproduced. 

We propose a scenario in which the cloud gravitationally 
contracts at significantly below the free-fall speed. A star clus- 
ter forms by fragmentation and accretion. The stars, especially 
the massive stars that dominate the luminosity, provide various 
feedback mechanisms to counteract the contraction. The addi- 
tional pressure, be it thermal, radiative, or turbulent, decelerates 
the contraction, then dominates first in the inner part, before fi- 
nally disrupting the whole cloud. SgrB2(M) is just beginning 
to drive out the gas, while large-scale global infall and probably 
also infall among very localized pathways in the interior remains 
ongoing. 

This study illustrates the importance of high-frequency lines 
in constraining the source structure, and demonstrates the great 
potential of HIFI, which delivers velocity-resolved spectra at 
these frequencies. Ground-based telescopes may not be able to 
reach high enough frequencies. In our case, detecting the HCN 
12-11 line was necessary in order to unambiguously trace the 
reversal of infall, so this result is unique to HIFI. We expect that 
a systematic study of high-J HCN transitions, even higher than 
the 12-11 line seen here, with a high signal-to-noise ratio, and 
toward a sample of hot core sources, would be highly rewarding 
in finding more sources in this particular stage of evolution. 
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